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We develop the statistical mechanics of the Hopfield model in a transverse field to investi- 
gate how quantum fluctuations affect the macroscopic behavior of neural networks. When the 
number of embedded patterns is finite, the Trotter decomposition reduces the problem to that 
of a random Ising model. It turns out that the effects of quantum fluctuations on macroscopic 
variables play the same roles as those of thermal fluctuations. For an extensive number of em- 
ly-j • bedded patterns, we apply the replica method to the Trotter-decomposed system. The result is 

summarized as a ground-state phase diagram drawn in terms of the number of patterns per site, 
\ a, and the strength of the transverse field, A. The phase diagram coincides very accurately 

with that of the conventional classical Hopfield model if we replace the temperature T in the 
' latter model by A. Quantum fluctuations are thus concluded to be quite similar to thermal 

fluctuations in determination of the macroscopic behavior of the present model. 

B 

g 

o 
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§1. Introduction 

X 

Statistical mechanics has been applied successfully to the analysis of various problems in neural 
networks. In particular, the Hopfield model, a prototype of associative memory, was solved explic- 
itly by combining well-established techniques in the mean-field theory of random spin systems. In 
the case of a finite number of embedded patterns^ the application of the mean-field method for 
the Ising ferromagnet to the neural network was shown to prove the existence of a memory-retrieval 
phase as well as mixed (confused) states. When the number of embedded patterns is extensive,!^ 
it is necessary to introduce the replica method in order to investigate the properties of the network 
at finite temperatures, similarly to the case of the Sherrington-Kirkpatrick model of spin glasses.!^ 
The resulting phase diagram is characterized by three macroscopic phases, namely, the retrieval, 
spin glass and paramagnetic phases. 
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The introduction of temperature in these studies was motivated by the apparent randomness 
in signal transmission at a synapseS A pulse reaching the terminal bulb of an axon does not 
always cause the release of neuro-transmitters contained in vesicles. The probability of release 
is generally small.!' & This randomness in signal transmission is usually taken into account as 
thermal fluctuations, which leads to the stochastic formulation of the problem in terms of the 
kinetic Ising model.& However, detailed considerations of the origin of this randomness suggest 
that quantum effects may be the major driving force to cause uncertainty in the release of neuro- 
transmitters from vesicles into the synaptic cleft. For example, StappS* pointed out that the 
migration of calcium ions in the bulb is a quantum diffusion process and thus the uncertainty in 
the ion positions leads to quantum fluctuations in the signal transmission at a synapse. Beck and 
EccleJ a rg ned that the nneertainty in the positions of hydrogen a,o m8 in the vesicnlar gr id is the 
origin of quantum fluctuations of comparable order of magnitude as thermal fluctuations in the 
brain at room temperature. These investigations indicate the necessity to treat randomness in the 
signal transmission in terms of quantum mechanics, not simply as thermal fluctuations as has been 
the case conventionally. 

The relation between quantum mechanics and the brain functioning has been discussed also in 
the context to clarify the fundamental significance of wave functions and observations in quantum 
mechanics (see Refs. [5] and |7] and references therein). However, few of the previous investigations 
in this area have paid attention to the behavior of macroscopic observables. Discussions have been 
given mostly in terms of microscopic wave functions, though experiments are often carried out on 
macroscopically observable variables. We should point out that the superposition of various micro- 
scopic state vectors does not always lead to the uncertainty in macroscopic observables. Quantum 
fluctuations may work in a manner similar to thermal fluctuations, leading only to weak stochastic 
deteriorations of observed values of macrovariables. We will show this effect explicitly in the present 
paper. 

We therefore have sufficient reasons to introduce quantum fluctuations into neural networks. It 
is in general difficult to reflect directly microscopic quantum processes in a simple model amenable 
to analytical investigations. We thus adopt the Hopfield model in a transverse field, 

w = -E J *x^-A5>r , (i.i) 

(u) « 

where af and of are the components of a Pauli matrix at site i. The interactions Jij are given by 
the Hebb rule as specified explicitly in §2. Admittedly this model is not a faithful reproduction 
of real processes in the brain. For example, the state of a neuron has quantum uncertainty in an 
eigenstate of the Hamiltonian (|1.1|), quite an improbable situation in reality. However, our purpose 
is not to explain the brain itself in detail. (One may argue in this regard that the Hopfield model 
without the transverse-field term is already inadequate as a model of the brain.) We rather aim 
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to clarify the role of quantum fluctuations in large-scale networks at a phenomenological level. We 
believe that the present system (1.1) serves as a first step toward this goal. 



This paper is organized as follows. In §2 the transverse-field Hopfield model (1.1) is solved for 



a finite number of embedded patterns. It is shown that quantum fluctuations have almost the 
same effects as thermal fluctuations on macroscopic properties of the network. The case with an 
extensive number of embedded patters is studied in §3. The resulting phase diagram turns out to 
be almost the same as that of the classical Hopfield model qualitatively and even quantitatively. 
(We call the model with A = in (1.1) the classical Hopfield model in this paper.) The last section 



is devoted to discussions. A preliminary report of a part of the present work has already been given 
elsewhere.!^ 

§2. Finite Number of Patterns Embedded 

We first consider the case in which the number of embedded patterns p remains finite in the 
thermodynamic limit. 

2. 1 Formulation 

The Hamiltonian of the Hopfield model in a transverse field has already been given in ( |1.1[) as 

W = - E J *°i a j - AX)o? = Wo + Wi. 

ft?) * 



The interactions Jjj obey the Hebb rule, 



1 p 



Jij = E tftf . ( 2 - 2 ) 
fj,=i 



with = 1 or — 1 randomly. The summation over the indices (ij) in (2.1) runs over all combinations 
of pairs of sites. The partition function of this quantum system can be represented in terms of 
simple Ising variables by the Trotter decomposition,!' 

f 

M— »oo V / iW-+oo 

where 



Z= lim Tr ( e -PH Q /M e -Wi/M\ M = Um z (23) 



(a M p M N \ 

~TT~fj E E E ZiCjViKVjK + B E °iK(Ti,K+l ■ (2.4) 
1V11 * K=l(ij)n=l K=\i=l j 

The coupling constant B in the Trotter direction is related to the coefficient A of the transverse 
field term in ( |2.lD 

B = -logcosec^^- • (2.5) 

2 M v ; 

We follow the standard procedure to decompose the double summation over (ij) using a Gaussian 
integral, 
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Z M = [ II dm ^M £ eX P ( "TiT £ + TF ££ ™>K^iOiK °iK<?i,K+l ) , (2.6) 

J Kfi o- \ Z1V1 Ky 1V1 Ky i Ki J 

where we have ignored the overall constant which is irrelevant for the following arguments. 

In the thermodynamic limit N — > oo with p kept finite, the saddle point of the integrand of ( |2.6D 
yields the equilibrium free energy per spin as 

/ = olJ £ m ^ " T « lQ g £ ex P 77 £ ^KvPok + BJ2 °k<tk+i » , (2.7) 

i<> a \ Kfj, K j 

where the double brackets <C • • • 3> denote the average over the randomness of embedded patterns 
We have assumed the self- averaging property of the free energy to derive the above expression. 
The saddle-point condition leads to the equation of state 

m K » =« » , (2.8) 

where the brackets (• • •) stand for the average by the weight 

\ Kn K j 

Equation ( |2.8D shows that the parameter represents the overlap of the spin configuration 

in the Kth. Trotter slice with the //th embedded pattern. The solutions of ( [2.8| ) describe the 
thermodynamic properties of the system. 

2. 2 Symmetric solution near the critical point 

Let us first discuss the symmetric solutions of the equation of state ( |2.8[) in the form = m 
for all K and \i < /, I being a given integer. If /i exceeds I, rriKa = 0. The stability of this type of 
solutions will be considered later. The free energy fl2.7|) is now written as 

/ = ^m 2 -T«log^exp(^^^^ + J BE^^+i) » ■ ( 2 - 9 ) 

cr V K £t=l K J 

To determine the critical temperature, we expand (|2.9[) to 0(m 4 ) as 



;oc 



/ = V -Tlo g Z - M « z f » { fe aK 



2 fa u 2M 2 

/3 3 m 4 ' x 



K 



«*f X £^ )oc , (2.10) 



24M 4 
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where 

a 

and 

In (2.10) the brackets (• • -)o c represent the cumulants calculated from the average (• • -)o defined by 
In order to evaluate the cumulants appearing in (2.10), we calculate the generating function 

z h ,M = J2 ex P \ B J2 aK(7K + 1 + ~nJ2 aK ) ■ ( 2 - 12 ) 

a \ K K ) 

In the limit of large M, we find, using the Trotter decomposition formula: 

Z h = lim Z hM = Tre haz+ P Aax 

M— >oo ' 



= 2cosh v //i 2 + /3 2 A 2 . (2.13) 
Expansion of this equation in powers of h gives 



log Z h = log Z 



h 2 tanh a 
+ 2^ 

M— >oo 



/z 4 (a — tanh a — a tanh 2 a) 
+ ~ IS? ' ■ (2 ' 14) 

with a = /3A. Comparison of (2.14) with (2.10) gives the cumulants in the limit M — > oo as 



1 ,/ V\ tanh a 
A^V^J )0c = ^' 

1 / \ 4 3 

■j^(y£2<7K J )oc = ^(a-tanha-atanh 



Using the relations 

< zf » = Z , (2.15) 

< z 4 » = Z(3Z - 2) , (2.16) 

we finally obtain the explicit expression of the free energy (2.10) as 

/ 1 2 (3m 2 tanh a 

— = -m • 

12 2 a 

, , , /3 3 m 4 a — tanh a — a tanh 2 a 

-(3/ -2)-^ -3 . (2.17) 
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The critical condition is thus expressed as 



tanh a 



or 



tanh/? c A = A 



(2.18) 



This formula shows that the critical temperature does not depend on the parameter I similarly to 
the case of the classical model.0) The critical line (2.18) is drawn in Fig. 1. The asymptotic form 



Fig. 1. The critical line separating the ordered and disordered phases. 



of the overlap order parameter around the critical temperature is derived from (2.17) as 



Pc( 3 ' - 2 )9{a c ) 



where o c = /3 C A and 



g(a) = — ^(a — tanh a — a tanh 2 a) . 



2.3 Symmetric solution in the ground state 
In the limit of large M, the free energy (|2.9|) is written as 



/ = hm 2 - T < log Trc&m'n+P*** > 

according to the Trotter decomposition formula. This expression is further evaluated in the limit 
T -» as 

/ = i/m 2 - T < log 2 cosh (3<J m 2 zf + A 2 > 



\lm 2 - < ^Jm 2 zf + A 2 > . (2.20) 
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The equation of state is obtained from (2.20) as 



//„ - i ; ^ » . (2.21) 

m 2 zf + A 2 



For a finite value of the overlap order parameter m, (2.21) reads 



I =< = > 

m 2 z 2 



i +A 2 

In consideration of (2.15), this equation implies that m decreases to as A approaches 1. Thus the 
critical value of A in the ground state is 1. The free energy (or the energy) at this critical point is 
— 1 independent of I as is apparent from (2.20). 

It is straightforward to derive the explicit forms of the order parameter and the free energy for 
small values of I from ( 2.21| ) and (2.20). The results are given by 

yi^A 2 ", / = _!(! + A 2 ) (2.22) 



m 
and 



2 



m = ivT^, / = -!(! + A) 2 (2.23) 



for I = 1 and I = 2, respectively, and 



3 1 
+ 



V9m 2 + A 2 Vm 2 + A 2 
/ = ^m 2 - I (V9m 2 + A 2 + 3y/m 2 + A 2 ) (2.24) 

for / = 3. These results are plotted in Fig. 2 for the order parameter and in Fig. 3 for the free energy. 
These figures show that in the ground state quantum effects represented by the parameter A play 
very similar roles to thermal fluctuations represented by the temperature T in the classical mo del) 
in which the order parameter and the free energy behave in almost the same way as functions of T 
as in Figs. 2 and 3 if we replace A by T. 

2.4 Stability of the symmetric solutions (I) 

It is necessary to check the stability of the symmetric solutions given in the previous subsection. 
The overlap order parameter rriKfj, should depend upon K and /i in general. The possible depen- 
dence of rriKfi on the Trotter number K is investigated in the present subsection. The pattern 
number dependence will be treated in the next subsection. 

Let us assume that m,Kn = ttik if fJ> < I and ra^ = otherwise. Then, the free energy fl2.7| ) is 
expanded to second order of mx as 



-TlogZo + -^^2m 2 K ~ ^2 ™>Km K >{o- K o-Ki)o = -TlogZ + A KK miKm K i 

K KK' KK' 



(2.25) 
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Fig. 2. The ground-state magnetization as a function of A for I = 1, 2 and 3. 



Fig. 3. The ground-state free energy (or the energy) as a function of A for I = 1,2 and 3. 



where 



A 



KK> 



$KK> 



M 



{vk&k')q ■ 



(2.26) 



The average (• • -)o was defined in (^7TT|) . The expansion of the free energy (|2.25[) corresponds to the 
approach to the critical point from the paramagnetic phase because the order parameter vanishes 
in equilibrium in the high-temperature paramagnetic phase. Correspondingly, from ( 2.26| ), all 
eigenvalues of the matrix A^k' are seen to be positive for small (3 or large T, implying the stability 
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of the paramagnetic solution mx = 0. The critical temperature is determined by the vanishing 
point of the minimum eigenvalue Ao of the matrix Arki, 

A = 1-4^ {<TK<TK')o 
KK' 



2 



tanh 3 A 

— ► 1 1— ■ 2.27 

AZ-oo A 

The critical condition tanh/? c A = A of ( |2,18| ) is thus recovered. An important point here is that 
the eigenstate corresponding to the minimum eigenvalue Ao is a uniform mode mx = Tn. This 
means that an ordered state uniform in the Trotter number K is formed slightly below the critical 
temperature. Therefore, we conclude that the symmetric solution rriKfi = m is stable against the 
Trotter-number dependence slightly below the critical temperature. 
At an arbitrary temperature the free energy is written as 

/ = ^E<- T « lo sE ex p(^E m ^K + <?K<JK + ^j » • (2.28) 

To investigate the stability of the symmetric solution, we check the eigenvalues of the Hessian 



d2 f 



drriKdrriL 



1 ' 6 kl-172 « z l ^kol) ~ (ctk)(vl)) » , (2.29) 



M M 2 



where the brackets (• • •) denote the average with respect to the weight 

6XP (fw~^ aK + B J2 a K a K+ij ■ 

Since the value of the matrix element Akl depends only on the difference \K — L\ and the quantity 
in the double brackets <C • • • 3> in ( |2.29| ) is positive, the lowest eigenvalue of this Hessian is given 
by 

M 
L=l 

The symmetric solution is stable if Ao is positive. Explicitly, 

Q 2 ( QziTFl \ 

MAo = /-/3^^<<logE-P^E^ + ^E^ + iJ>> 
B 2 

— ► I - T-^r < log Tr e P z i m °z+P A(T x > 
M^oo dm 2 



I - A 2 < - 2 1 tanhdJm 2 z 2 + A 2 > 

(m 2 zf + A 2 ) 3 / 2 v ' 
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r 



pm 2 <C 9 ' A9 sech 2 /?Jm 2 z, 2 + A 2 > . (2.30) 
m 2 zf + A 2 v 1 



At T = 0, the above expression reduces to 



.2 



MX = l- A 2 « Zl » 

(m/zj + A 2 ) 3 / 2 



For A — > 0, MXq tends to I, a positive value. If, on the other hand, A is close to 1 (A = 1 — e 
with e <C 1), MXq approaches 2el, again positive. We thus expect that the eigenvalue Ao remains 
positive between these two limiting values of A when T = 0. 

In the case of general finite temperature, we have numerically confirmed the positivity of Ao for 
1 = 3. It is expected that the same property holds for other values of I, though it is difficult to 
prove it explicitly. 

2.5 Stability of the symmetric solution (II) 

We next check the pattern number Qu) dependence of the solution of the equation of state. If we 
assume m^ = m M for /j, < I and m/^ = otherwise, the free energy ( |2.7| ) reads 

/ = ^E m ^- r<<1 °gE ex p ( TiJ2 m ^J2 a K +B J2 a ^K+i) » 

M \ 1V1 H K K J 

— ► -m 2 - T « log 2 cosh pJ (m • £) 2 + A 2 > . (2.31) 

M^oo 2 v 



The vector m has components (mi, m%, • • • , m/, 0, 0, • • •), and similarly for £. The Hessian around 
the symmetric solution is given by 



d2 f 

A 



dmndrriy 



m a =m 



^ " A <K (m 2 z 2 + A 2 )3/ 2 tanh N m2z ? + A2 » 
m 2 z, 2 + A 2 " 



-/3m 2 <C ff^r 9 sech Wm 2 z 2 + A 2 > . (2.32) 
m 2 zf + A 2 v ' 



The diagonal element of {A^j,} is written in the following form 



tanh /3y m 2 zf + A 2 
(m 2 zf + A 2 ) 3 / 2 



ei = A w = 1 - A z < - - v 9 : a9xq/9 » 



- /3m 2 « ^ ^ » , (2.33) 
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where 



m 2 z 2 



Q = < o 9 ' ao tanh2 P\ m 2 z? + A 2 > 
ra^zf + A 2 v « 

= «(^> 2 > . (2.34) 
All off-diagonal elements have the same value: 

e 2 = \u ((* ¥= v) 

= " A2 « {m 2 z ff A 2 )3/ 2 tanh/J^^ + A 2 » 

_/?m2 « m 2 1f+A 2 SeCh2/? V /m2 ^ + A2 » • ( 2 - 35 ) 

From (2.32), (2.33) and (2.35), it is easy to see that there are three eigenvalues of the Hessian 
matrix {A^}, 

Ai = e x + (Z - l)e 2 , (2.36a) 
A 2 = e 2 , (2.36b) 
A 3 = ei - e 2 , (2.36c) 

with the degeneracy p — 1,1 and Z — 1, respectively. 

To investigate the stability of the symmetric solution near the critical temperature, we expand 
various terms in (2.33) and (2.35) to second order of m using (2.34): 

9 z? m 2 l 

171 < 2 2 7 A2 »~ "a" ' 

m z zf + A z A 

2 , tanh 2 B A 
/ _ 

and 



tanh/?wm 2 2: 2 + A 2 \ ?>m 2 l 8m 2 ] \ 

« 7 2 2 A L/ 2 »« ( tanh /? A " St tanh ^ + OA , 2/3A • (2-37) 
(m 2 z 2 + A 2 ) 3 / 2 A 3 ' 2A 2 2Acosh 2 /3A/ 



These equations and ( |2.19| ) yield 



2 /tanh/3A A 2(1 - A 2 ) 

61 ~ 37^2 [ A J ~ 3Z-2 6 ' 

6(A" 1 tanh/3A- 1) 3(1 - A 2 ) 

62 ~ 37^2 ~ 3Z-2 6 ' 

where e represents the deviation of the temperature from the critical value, 

e = T C (A) - T . 
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The signs of the eigenvalues in (2.36) are then found as Ai, A2 > and A 3 < 0. Since the degeneracy 
of the last eigenvalue A3 is Z — 1, the symmetric solution with I > 2 is unstable against fluctuations 
in the direction of the eigenvector corresponding to A3 slightly below the critical temperature. This 
property is exactly the same as in the case of the classical Hopfield model.0) 
When T = 0, the matrix elements (2.33) and (2.35) reduce to 



ei 



A 2 <C 



-A 2 



(m 2 zf + A 2 ) 3 / 2 



(m?zf + A 2 )3/2 

For I = 1, Ai = e\ = 1 — A 2 > 0, implying the stability as expected trivially for the one-component 
solution. When I = 2, the eigenvalues are given by 



Ai 
A 2 
A3 



A 2 

2 
1 

A 



1 

2A 



(2.38a) 
(2.38b) 

(2.38c) 



The first eigenvalue Ai is always positive in the range < A < 1, while the second one A2 is positive 
for A > 0.618 and is negative for A < 0.618. The third eigenvalue A3 is always negative. Thus, 
this I = 2 solution is unstable at T = for any A between and 1. 

The symmetric solution with I = 3 has the following eigenvalues for T = 0: 



Ai 
A 2 
A3 



1 



A'- 



+ 



1 



A 2 
1 

4 



1 



( 9m 2 + A 2)3/2 ( m 2 + A 2)3/2 

1 3 



+ 



(9m 2 + A 2 ) 3 / 2 (m 2 + A 2 ) 3 / 2 
A 2 



(m 2 + A 2 ) 3 / 2 ' 

These eigenvalues are plotted as functions of A in Fig. 4. The first and second eigenvalues Ai and A2 
are both positive, but A3 changes its sign at A = 0.494. Similarly, the behavior of the eigenvalues 
for I = 4 and I = 5 is shown in Figs. 5 and 6. This analysis for I up to 5 thus suggests that the 
even-/ solution is always unstable while the odd-Z solution is stable in a finite range < A < A c 
when T = 0. 

To confirm this conjecture, we have expanded the eigenvalues Ai,A2 and A3 as functions of 
e = 1 — A for small e. The equation of state ( |2.21 ) has the solution m = c^ft with c = y/2/(3l — 1). 
This gives 

1 



ei 

C J 2 



:(3lc 2 -2)e , 



3c 2 e . 
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Fig. 4. Three eigenvalues of the Hessian for I = 3 and T = 0. The third eigenvalue A3 crosses at A = 0.494. 



Fig. 5. Three Hessian eigenvalues for I = 4 and T = 0. 



Then, the eigenvalues (2.36) are 



Ai = 2e , 
As 



3Z - 2 ' 
3/ + 2 



3/ - 2 

Since the third eigenvalue A3 has degeneracy I — 1 and is negative, we conclude that all symmetric 
solutions with I > 2 are unstable near A = 1. This, together with the analysis of eigenvalues for 
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Fig. 6. Three eigenvalues of the Hessian for I = 5 and T = 0. The third eigenvalue A3 crosses at A = 0.426. 



I = 2 to 5 explained before, suggests that the coefficient A of the transverse-field term of the present 
model in the ground state has effects very similar to the temperature of the classical Hopfield model 
(A = 0) in which the even-/ solutions are unstable for any T while the odd-/ solutions are stable 
below certain temperatures. & 

The instability of the even-/ solutions suggests the existence of asymmetric solutions of the 
equation of state ( |2.8| ) as was the case in the classical model-S-* We actually have found many 
asymmetric solutions. An example of a solution of the type (m, m,u, u, u, 0, 0, • • •) is shown in 
Fig. 7. Although we were not able to find general rules for the existence and behavior of asymmetric 
solutions, we observed that, when T = 0, A has effects very similar to T in the classical model. 

§3. Extensive Number of Patterns Embedded 

In this section the number of patterns embedded, p, is assumed to be proportional to the system 
size N. We closely follow Chap. 10 of Ref. [l^ in our presentation. 

3. 1 General form of the free energy 

According to ( |2.6| ), the replicated partition function averaged over the quenched randomness is 
written for a finite Trotter number M as 



N(3 
2M 



E 

Knp 



m Kfj,p 



+^ J2 J2 m KfMp ^a ipK +BJ2 <?ipKcri p ,K+i ] » , (3-1) 



Kp,p i 



Kip 



14 



Fig. 7. An asymmetric solution of the type (m, m, u, u, u, 0, 0, • ■ ■) for T = 0. 



where p (= 1, ■ • • , n) represents the replica index. The overall constant is irrelevant and is ignored 
here. 

We consider the situation in which a finite number of patterns have nonvanishinig overlaps. That 
is, rriKup is of order unity for /j, = 1, 2, • • • , s and is of order 1/y/N for other ^'s with s being a 
finite number. Then, the configurational average <C ■ ■ ■ S> in (3.1) can be evaluated explicitly for 
patterns fx > s. For the fith pattern with fx > s, 




(3.2) 



The terms containing {rriKfip}p>s in (3.1) and (3.2) are collected in a single formula, 




where 



A-Kp,Lcr = $Kp,La 



NM 



^VipK&iaL ■ 



(3.3) 



The integration of over {rriKpp} gives 



[J drriKppEfj, = const, x (det A) 

Kp 
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The product of this result over p = s + 1, ■ ■ ■ ,p is 

(detA)-( p - s )/ 2 « (detA)- p / 2 

V 1 Kp 

Here we have neglected the finite number s in the exponent since it is small compared to the 
extensive number p. The eigenvalues of the matrix A are denoted as \k p - 

The eigenvalue Xk p depends on the spin configuration as is apparent from ([T^). To avoid this 
complication, we introduce another matrix A without explicit dependence upon the spin configu- 
rations: 

P P 
^K P ,La = S Kp ,La - —q pa (KL) - 5 pa —S p (KL) . (3.4) 

This matrix A is equal to A, if 

q pa (KL) = { N i (3.5) 
(p = a) 

S P (KL) = j-fYl a ipK<7ipL (3.6) 

i 

are satisfied. Physically, q pa (KL) is the spin glass order parameter for spins at the Trotter slices 
K and L, and S P (KL) is a measure of quantum fluctuations. If there are no quantum fluctuations, 
the spin configuration Oi P K does not depend on the Trotter number K and S P (KL) = 1. Quantum 
fluctuations reduce S P {KL) from unity. Hence 1 — S P (KL) gives the scale of quantum fluctuations. 
It is convenient to rewrite a function of \k p as 



and 



G{\ KP }= I n dq pa (KL)n n &s p (kl) 

(Kp,La) P (KL) 

xS ^q pa (KL) - ^^aipKcricrLj 



1 



x5 ( S P (KL) - — dpKOipL J G{X Kp } , (3.7) 



where (KL) denotes an arbitrary combination of K and L including K = L, and (Kp, La) expresses 
an arbitrary pair except for p = a. Using the Fourier representations of the delta functions, (3.7) 
is written as 



G{X Kp }= f II dq pa (KL)dr pa (KL)l[ ]J dS p (KL)dt p (KL) 
(Kp,La) P [KL) 



x exp 



Na(3 2 _ r , , TjrT . a(3 2 , T ^ T ,^ 
■—j^-r pa (KL)q pa (KL) + -^r pu (KL) ^ (T ipK (T iuL 
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A ( P \(KL)S P (KL) + j^t p (KL) ]T a ipK a ipL 



M 2 



G{\k p } , (3. 



where we have ignored the overall constant. 
The average of the replicated Z is thus 



«Z"»= f l[dm Kpp [] dq pa (KL)dr pa (KL)l[ [] dS p {KL)dt p {KL) 
Kpp (Kp,La) f (KL) 

N 
~2 



x exp 



■ 2M 2^ L m K W ~ — E lo § A ^ 



iVa/3 2 

-jj^r -pffv— i 2M 2 

(K p,Lcr ) P (A-L) 



E r pa (KL)q pa (KL) 



Nap 2 



E E t p (KL)S p (KL) 



x < exp 

<7 



M 



E E m *>P E C°V^ + 5 E E a ipKVip,K+l 
i i Kp 



Kp p,<s 
(2 

2Af2 



+ ^2 2 E r pAKL)a ipK a iaL + |^ EE E tp(KL)a ipK a ipL 



i (Kp,La) 



i P (KL) 



> (3.9) 



If we write the above integrand as exp(— iV/3/), the saddle point of the integral yields the free 
energy per spin / as 



1 



Kp p<s ' Kp 

(Kp,La) P (KL) 



- T < logE ex P 



M 



E E m Kpp^<? pK + £ E a pK<*p,K+\ 



Kp p,<s 



2M 2 



(Kp,La) 



Kp 

+ tttL E r pc{KL)a p K(T aL + E E t p (KL)a pK a pL 



P (KL) 



> . (3.10) 



We have used the self- averaging property of the free energy in dropping the site index i in the above 
equation. The stationarity conditions of the free energy give the equations of state: 



mKpp =< ^{<? p k) > , 
q pa {KL) =< (a pK )(a aL ) > , 
r pa (KL) = - E < m K^ P m Lp(T > , 



(3.11) 
(3.12) 
(3.13) 
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S p {KL) =« (a pK a pL ) » , (3.14) 
t p (KL) = - ^2 < m Kpp m Lpp > . (3.15) 

/i>S 

The physical meaning of the parameters m^fip, Qpa(KL) and r pa (KL) is the same as in the classical 
Hopfield model: mxpp is the overlap, q pa (KL) denotes the spin glass order parameter and r pa (KL) 
represents the effects of uncondensed patterns. The deviation of S P (KL) from unity reflects quan- 
tum fluctuations as mentioned before. The last quantity t p (KL) is for effects of uncondensed 
patterns in the same replica, or in other words, the diagonal element of r pa (KL). 

3.2 RS solution in the static approximation 

It is very difficult to solve the equations of state in their general forms (3.11)-(3.15). We instead 
look for the solutions in the replica symmetric (RS) subspace under the static approximation.00) 
We thus neglect the dependence of the order parameters on the replica index (the RS approxima- 
tion) and the Trotter number (the static approximation): 

m K pp = m p , (3.16) 

q pa (KL) = q , (3.17) 

r pa (KL) = r , (3.18) 

t p {KL) = t , (3.19) 

The stability of the replica symmetry will be considered later. The static approximation is expected 
to give at least qualitatively reliable results as long as the parameter A representing quantum effects 
is not too large as was the case of the SK model in a transverse neldSIll) The consistency of the 
RS and static approximations can be checked from the view point of another approximate method 
as explained in §3.3 and Appendix A. 
The free energy (3.10) divided by n (the number of replicas) is now 
1 9 a v-^, , ad _ ,o. „. aB _ o _ 

' Kp 



T 

■— <C log exp 
n _ 



(3 

Kp fi<s Kp 

EE^) 2 -e(e^) 2 } + ^e(e^) 2 » -(3-2D 

The summation of logA^p appearing above is carried out as follows. There are three values of the 
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matrix elements of {^-Kp,La} denned in (3.4), 



-■jLs (p = a,K^L), 
( P = <t,K = L), 

under the static approximation. The eigenvalues of this matrix are easily found to be 

M - 1 



Ai = 1 - — 

M M 

13 M-l 
A -1 + P 9 



0S - (n - l)(3q , 



with degeneracies l,n — 1 and (M — l)n, respectively. Thus, we have 



lim - V log X Kf 

n-»0 n„ 
Kp 



i _ J_ _ M-l 

M M 



/35 + /3g 



+ (M-l)log(l-A + A 5 ) 



0q 



log(l-/35 + /?g)- - -P + 0S. 
M^oo 1 — pb + pq 



(3.22) 



The Gaussian integral enables us to decompose the last term in the exponential appearing in (3.21) 
into independent replicas: 



exp 



ap 2 
2M 2 

Dz exp 



mEE^ + (*-»•) E E°>* 

V p A' / p V A y 



M 



E°>* + 



pA 



a(3 2 (t - r) 
2M 2 



E E 



p \ A 



(3.23) 



where Dz denotes the Gaussian measure e 2;2 / 2 dz/v / 2vr. Then, the summation over a p in (3.21) 
can be carried out independently for each p as 



log E ex P 



M 



EE^E "v^ + 5 EE <v^p,#+i 



+^ r (EE^l + ^(*-^e(e^1 



P A 



p \ A 
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= log J Dz j^exp 



JL 

M 



K 



K 



K 



K 



i Q n J D-g log exp 



K 



K 



K 



K 



n 



J Dzlog^] J Dw exp 



K 



K 



K K 



P 



M^oo 



n 



J Dz log J Dw Tr exp 
J Dzlog J Dw2cosh/3W 



m ■ £ + y/arz + y a(t — r)io 



1 2 



+ A 2 



(3.24) 



The total free energy is then obtained from (3.21), (3.22) and (3.24) as 



, 1 2 a 



log(l -PS + /3q) 



Pq 



l-PS + Pq 



- P(i - S) 



+ f(tS-rq) 



-T <C J Dz log J Dw 2 cosh pJ m • £ + \/arz + \j a(t — r)w 



+ A 2 > . (3.25) 



3. 3 Equations of state at T = 

Variation of the free energy (3.25) gives the equations of state for the order parameters. Let 
us consider the case in which only one of the patterns is retrieved, m M = S^i ■ m. For brevity of 
expressions, we introduce the following notations, 

g = m + \farz + \j a(t — r)w , 
u = ^Jg 2 + A 2 , 
Y = J Dw cosh Pu . 

The equations of state obtained by the variation of the free energy (3.25) with respect to m, q, s, r 
and t are, respectively, 

m = JdzY^J Dwgu^smhPu , (3.26) 
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1 -/»+*• 2 <3 - 28) 

q = J Bzfy- 1 J Bwgu' 1 sinh pv}j , (3.29) 
< S = y Dzy^ 1 ^ Du> g 2 u~ 2 cosh /fa 

+TA 2 / Bw u~ 3 sinh /fa j . (3.30) 



Quantum effects play most important roles in the ground state where thermal fluctuations are 
absent. We therefore only consider the case T = hereafter. The equations of state simplify 
considerably in the T = limit because q = S and t = r when T = as proved below. For very 
small T, we find from (3.29) and (3.30) 



S = J BzY- 1 J Bw g 2 u- 2 cosh (3u 

> J BzY' 1 J Bw g 2 u~ 2 sinh /fa 

> J Bz (y' 1 J Bw gu^ 1 sinh (3v)j = q . 



If we assume that S is strictly larger than q at T = 0, (3.28) immediately leads to t = r and also 
(3.27) to r = because —f3S + (3q diverges. Then, u is a constant (m 2 + A 2 ) 1 / 2 , and hence (3.29) 
and (3.30) imply S = q, a contradiction to the assumption S > q. Therefore, the equality S = q 
holds in the ground state. 

To derive the ground-state equation of state for r from (3.27), we have to estimate the T — > 
limit of 0(S - q). Comparison of (3.29) and (3.30) leads to 

lim (3{S -q) = A 2 f Bz- — — ee C . (3.31) 

J [{m + ^Fz) 2 + A 2 ] 

The equations of state at T = are thus written as 

Vz -. V , (3.32) 

. i r7ZZ-,\ 2 i A 2 



\J {m + \farz) 2 + A 2 

/ / \2 

(m + varz) 
(m + y'aFz) 2 + A 2 



(3.33) 
(3.34) 



(i-cr ■ 

The same equations of state can also be derived by a direct mean-field analysis as explained in 
detail in Appendix A. The reason why we have presented the replica analysis in the present section 
is three fold. First, the replica method gives the free energy by means of which we can distinguish 
two different retrieval phases as explained in the next subsection. Second, the AT line, namely the 
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stability limit of the RS solution, can be determined by the replica formalism as discussed in §3.5. 
Lastly, it is useful to confirm that two different methods lead to the same results so that we acquire 
confidence in the appropriateness of the present approximations. 

3.4 Phase diagram at T = 

To draw the phase diagram, the equations of state (3.32), (3.33) and (3.34) have been solved 
numerically. The result is drawn in Fig. 8. The asymptotic form of the overlap order parameter 



Fig. 8. The ground-state phase diagram of the Hopfield model in a transverse field. "R-I" stands for the retrieval 
phase in which the retrieval states are the global minima, and "R-II" denotes the retrieval phase where the retrieval 
states are the local minima. The dashed line represents the AT line. 



around the The three phases are characterized by the relations m = q = (paramagnetic, P), 
m = 0, q > (spin glass, SG) and m ^ 0,q > (retrieval, R-I and R-II), respectively. The retrieval 
phase is separated into two parts, one with /r < /sg (R-I) and the other with /sg < /r (R-II), 
where /r is the free energy of the retrieval state and /sg is that of the spin glass state. 

The transition between the paramagnetic and the spin glass phases is of second order, and so the 
shape of the boundary can be determined analytically by expansion of the equation of state (3.33). 
The result is 

A = 1 + \[a . (3.35) 

This is exactly the same relation as the corresponding classical phase boundary if we replace A 
by T.B The other phase transitions (between SG and R-II and between R-I and R-II) are both of 
first order. Therefore, it is in general impossible to obtain the analytic expressions of these phase 
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boundaries. However, when a is very small, we can determine the asymptotic forms of the phase 
boundaries by expansions of the equations of state and the free energy as in the case of the classical 
model.!-* Details are described in Appendix B. The result is 



A ~ 1 - 1.95</a (3-36) 



for the transition between the SG and R-II phases, and 



33 _ 

A~l ^= 1 - 2.0625^/0 (3.37) 
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between the R-I and R-II regions. The relation ( 3.36j ) exactly agrees with that of the classical case 



if we replace A by T.EP The relation ( 3.3?1 ) shows that the R-I region of this model in the vicinity 



of (a = 0, A = 1) is wider than that of the classical case, T ~ 1 — 2.6-^/a,!^ though the deviation 
from unity is also proportional to y/a. 

When A = 0, the R-II phase changes into the SG phase at a = 0.1379 as it should @* The 
boundary between the SG and R-II phases is slightly reentrant in the low-temperature region. 
That is, for a fixed a slightly larger than 0.1379, the SG phase once changes to the R-II phase as 
A is decreased but the SG phase once again becomes stable for very small A. The same is true for 
the boundary between the R-I and R-II regions. The reentrance in these two cases is observed also 
in the classical model. It should be noted, however, that the effect of replica symmetry breaking 
treated in the next subsection obscures the significance of reentrance within the RS solution near 
a = 0.1379, again the same situation as in the classical model© 

3.5 AT line 

The stability of the replica symmetric solution against replica symmetry breaking can be checked 
following the standard procedure.!' 0* We keep the static approximation intact and investigate only 
the effects of replica symmetry breaking. Calculations are somewhat involved but straightforward. 
Details are given in Appendix |~C] . The stability limit of the replica symmetric solution (the AT 
line) is found to be given by 

a = arA / Dz 



(m + \/afz) +A 2 . (3.38) 

This AT line is also drawn in the phase diagram in Fig. 8. 

The region near a = 0.1379 and A = is drawn enlarged in Fig. 9. The AT line merges with 
the boundary between the R-I and SG phases at a ~ 0.1378 and A ~ 0.022. Since the replica- 
symmetric solution is unstable below this line, the reentrant behavior mentioned in the previous 
subsection may be an artifact of the replica symmetric solution. 

§4. Summary and Discussion 

The Hopfield model in a transverse field has been introduced and solved. For a finite number of 
embedded patterns, we have used the method of the Trotter decomposition to reduce the quantum 
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Fig. 9. The phase diagram near a — 0.1379 and A = 0. 



problem to a classical form. We found that the solutions of the equations of state in the ground 
state have quite similar properties to those of the classical model at finite temperatures. 

If the number of embedded patterns is proportional to the system size, it is necessary to employ 
the replica technique, in addition to the Trotter decomposition, to trace out the quenched random- 
ness. We cannot obtain the full exact solution and have to appeal to the replica symmetric (RS) and 
static approximations. The stability analysis of the RS solution leads to the AT line, above which 
the RS solution is at least locally stable. On the other hand, the static approximation probably 
does not give the exact solutions for all values of the parameters except for the limits of A — > 
(the classical model) and a — > (the finite-p model). Nevertheless, the experience in the transverse 
SK model© suggests that the static approximation is expected to capture the qualitative features 
of spin systems in transverse fields. The resulting phase diagram has three phases, the retrieval, 
spin glass and the paramagnetic phases. The shapes of phase boundaries turn out to be almost the 
same as those of the classical model if we replace A by T. 

We have found that the quantum fluctuations have almost the same effects as thermal fluctuations 
in the Hopfield model. The uncertainties in signal transmission at a synapse have conventionally 
been treated in model analyses in terms of thermal fluctuations. The results of the present study 
show that we may instead consider quantum fluctuations without changing conclusions on the 
macroscopic behavior of the network. 

The above-mentioned fact means that quantum uncertainties do not lead to many possible alter- 
native values of macroscopic variables but rather cause only quantitative deteriorations as thermal 
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fluctuations do. Quantum uncertainties at microscopic levels are far from the source of the simul- 
taneous existence of macroscopically distinguishable states. This point seems to be disregarded in 
some of the arguments concerning the significance of quantum mechanics in the functioning of the 
brainll 
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Appendix A: Equation of State by the Mean-Field Theory 

In this Appendix we derive the equations of state (3.32), (3.33) and (3.34) by a direct mean-field 
analysis.© The Hamiltonian of the Hopfield model (|2.1| ) with the Hebb rule (|2.2| ) is approximated 
by the following single-site effective Hamiltonian: 



Wi = -0f£3H-Aaf , (Al) 



where stands for the overlap 



m » N 



^ £#<«#• ( A - 2 ) 



The thermal expectation value of er| appearing in (A2) can be calculated easily under the Hamil- 
tonian ( |Aip to give 



U7 + A 2 



(af) = , / tanh/Vu 2 + A 2 , (A-3) 

where 



^ = £m^f. (A4) 

Let us now assume that a single pattern, the first one (fi = 1) for instance, is retrieved. Then, 
mi = m is of order unity and all other m^'s are of order 1/yN. The equation of state ( |A2| ) for 
[i = 1 is then written as 

1 » - m + yfarzi 



rn 



£ 



■ (m + \farzi) 2 + A 2 



x tanh/?y (m + \/arzi) 2 + A 2 , (A-5) 

where 

ii>2 

Since z% is proportional to the sum of extensively many random variables, it is not unreasonable to 
assume that z% is a Gaussian random variable according to the central limit theorem. Strictly speak- 
ing, the terms appearing in the definition ( |A6 ) are mutually correlated through the ^-dependence 
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of and the central limit theorem is not applicable. Nevertheless, the results obtained under 
this approximation of Gaussian distribution agree with those from the replica method in the case 
of the classical Hopfield model© We thus adopt the same approach in the present quantum case. 

The average value of z% is expected to be vanishing because £f is ±1 with equal probability. 
The variance of z% is unity if we define r in the following way: 

« E m ^i$ »=« E m l »= ar ■ ( A - 7 ) 

\Ai>2 / m>2 

Then, in the thermodynamic limit, the summation over i in (A-5) reduces to the average over the 
Gaussian distribution as 



, „ m + \/arz 

m = uz- 



(m + y/arz) 2 + A 2 



x tanh/3y (m + \farz) 2 + A 2 , (A-8) 

where T)z denotes the Gaussian measure e~ z2 ^ 2 dz/y/27r. 

To derive the equation of state for r, we need to evaluate m„ for v > 2 according to ( |A-7| ). For 
this purpose we decompose Uj defined in ( |A-4| ) into three parts by extracting contributions from 
// = 1 and n = v: 

Uii] = m + yfarzi + rf(m v , (A-9) 

with r]% = Although the definition of Z{ here does not include the contribution from \x = v in 

contrast to (A-6), this difference does not affect the results because the only property of Zj we use 



is its Gaussian distribution, which remains the same if we drop a single term from (A-6). It should 



be noted here that m u is of order 1/yJV and is quite small compared to the other two terms in 

(ED- 

We insert ( |A-9|) into the right hand side of (|A 2| ) with ( A-3 ) taken into account and expand the 
result to first order of m„ to find 



m " = ~\r E / tanh/^f 2 + A 2 



N 

N ^ 



+ A 2 



, 2 A A2N , /2 tanh/^f + A 2 + /\ 2 sechW?; 2 + A 2 
(uf + A 2 ) 3 / 2 v 1 vf + A 2 v 1 



(A-10) 



with Vi = r\\ (m + y/arzi). We now define q and 5 as 

9=^E^^^nh 2 /V^7^ 



?; 2 



Dz^ -r tanh 2 /? vV + A 2 (A-ll) 

w 2 + A 2 
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TA 2 



„ . - tanh/3\/ v 2 + A 2 
(t> 2 + A 2 ) 3 / 2 

with v = m + ^Jarz. Equation (A- 10) can then be solved for m u as 



(A-12) 



1 



l-pS + /3q 

1 ^ Vi 



E 



iV V v? + A 2 



tanh/3Jw 2 + A 2 . 



(A-13) 



From (A- 10) and (A- 11), we find 



(1-PS + (3qf < m l » 
p — 1 



!^>2 



E 



tanh 2 /?,/t; 2 + A 2 



^2 ^„2 + A 2 

= ag . 

Equation ( |A-7|) and the above relation lead to 

It is straightforward to check that the zero- temperature limits of (A-8), (A-ll), (A-12) and 
agree with ( |3.31| ), (3.32), (3.33) and (3.34) derived by the replica method, respectively. 

Appendix B: Phase Boundaries Near a = 

In this Appendix we derive the asymptotic forms of phase boundaries around the point {a = 
0, A = 1). The first boundary is between the SG and R-II phases and the second is between the 
R-I and R-II phases. 

To derive the boundary between the SG and R-II phases, we expand the right-hand side of the 
equation of state (3.32) assuming m and r are small: 



m = — 
A 



1 



m 



2A 2 



(m 3 + 3mar) 



Let us write 1 — A = e. Then, the above equation is approximated to order e as 

m + 3ar = 2e . 



(B-l) 



Similarly, the equations of state for C, q and r (( |3.31[ ), (3.33) and (3.34)) behave asymptotically 

as 

3 



C=(l + e) 

'2 



1 (m + ar) 

2 V ' 



q = m + ar , 

4(m 2 + ar) 



r 



(2e — 3m 2 — 3ar) 2 



(B-2) 
(B-3) 
(B-4) 



27 



To erase r from (B-l) and (B-4), we define x = m 2 + ar and write the above equations as 

X = m + (3x-2e) 2 ' (B - 5) 
2 

x = _( e + m 2 ) . (B-6) 
3 



Prom these equations we obtain 



2 / 2x 2 2a(e + m 2 ) 

-(e + m 2 ) =m 2 + (B-7) 
Using new parameters y and r defined by 

y = m 2 /^ , (B-8) 

r = e/x/a 7 , (B-9) 

( p-7| ) is expressed as 

~y 3 - T y 2 + y + T = , (B-10) 

which coincides with (5.12) of Ref. 0, though the definition of y is not the same. Since bifurcation 
occurs on the SG-R-II phase boundary, the derivative of this equation is also vanishing on this 
boundary, that is, 

^y 2 -2ry + l = . (B-ll) 
Eliminating r from (B-1C) and ( B-ll| ), we have 

b 3 + y = |y 2 + i 
y 2 - 1 2y ' 

and therefore 

2 5 + ^ 
7/ = and t = 1.95 • • • , 

or equivalently, 

A ~ 1 - 1.95ya . 

The boundary between the R-I and R-II phases is the line on which the free energy of the solution 
with m ^ and that of the solution with m = are equal to each other. We first derive the explicit 
form of the free energy at T = 0. All we have to do is to take the limit j3 — * oo in the expression 
of the free energy at finite temperatures (3.25). However, we should be very careful in taking this 



limit because both (S — q) and (t — r) are of order T, as shown in (3.28) and ( |3.31| ). In order to 
evaluate the last term of (3.25), — T <C • • • 3>, we expand this term up to first order of (t — r) for 
the single retrieval case to obtain 



-T < • • • > 



28 



7 Bziog I 



Dw 2 cosh 



„ r, i — ,2 77 / , (to + Jarz) Ja(t — r)w 

\ (to + yarzj + A z 



T{ Dz loo 



exp ( /3y (m + y / arz) + A 2 J / Dw exp /3 



, (to + y/orz) >/a(t — r)w 



m + \farz) + A 2 



+ exp I — (3y (to + y/arz) + A 2 I / Dw exp — /3 



m + y / arz) + A 2 
(to + y/arz) 2 



— T / Dz log 2 cosh /3 v (to + \farz ) +A 2 +-/3 2 a(t-r) / D; 

[/ V VV ; / 2 M ; 7 (to + a/otz) + A 2 

-» - / Dz-y/ (m + v^-z) 2 + A 2 - -a/3(i - r)g . 



(B-12) 



Substituting this formula to (3.25) and taking the limit T — > 0, we have 

H-" 2 -f( i+ T^- rc 



J ~Dz\J (to + \faxz) 2 + A 2 . 



(B-13) 



Next, we expand this expression with respect to to 2 and e = 1 — A. Note that the last term of 
(B-13) is not multiplied by a smallness parameter a/2 in contrast to the preceding term, and thus 
we have to calculate up to second order of e in this term. We aim to express the free energy only 
using the parameters y and r as ( B-10|) . For this purpose, we write all the parameters in terms of 
to 2 and e. In the solution with to 7^ 0, QB-1|) , (B-3) and (B-2) reduce to 

2 



i 2 

2 1 



q = -(to + e) 



ar 
C = 



1 + e 



-TO 



ttT + e 



(B-14) 

(B-15) 
(B-16) 

II). Substituting the 



On the other hand, in the solution with to = 0, we use (3.34) instead of 
relation q = ar to (3.34), we have (1 — C) 2 = a. Because the term 

exists in the original expression of the free energy (3.25), the condition C < 1 should be satisfied 
and we obtain 

2 



q = ar 



+ Va) 



C = 1- V« 



Using (B-8), (B-9) and (B-13)-(B-18), after long but straightforward calculations we have 



a 



/(to) - /(0) = - 



-v 

T 



2r [y + - + 4r + 1 



. 



(B-17) 
(B-18) 

(B-19) 
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Solving this equation together with ( |B-10| ), we obtain 



or equivalently, 
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y = 3 and r = — = 2.0625 
y 16 



A ~ 1 - 2.0625\/a . 



Appendix C: AT Line 

We derive the expression of the AT line ( |3.38| ) at T = 0. We follow the Appendix B of Ref.|| in 
the first half of the argument. To obtain the second derivatives of the free energy by the variables 
q pa and r pa under the static approximation, we retain only the relevant part in (3.10): 



f = jp Yl lo s x k p + r p°yp° 



Kp 

T -C log 2J exp 



pa 



M 



J2 £ <vr 



PA' 



K 



a/3 2 



+B J2 <r P K<T p ,K+i + ^rj2 H o P K<y aL 

Kp pa KL 



» . 



(C-l) 



We have set r pp = t„ in anticipation of the ground state relation r = t at the replica-symmetric 
point as explained in §3.3. 
The Hessian, or the second-derivative matrix of the free energy, is written as 



where 



d 2 f 



d 2 f 

dr a/3 dr 7 s 
d 2 f 



(C-2) 

(C-3) 
(C-4) 
(C-5) 



There is no difficulty in evaluating the derivatives written above by taking into account that {Xk p } 
is a set of eigenvalues of the Mn x Mn matrix ( [3.4| ) under the static approximation. Note that the 
replica-index dependence of the matrix elements should be kept untouched until differentiations 
are carried out. The result is almost the same as that of the classical Hopfield mo del§ For the 
submatrix A, 



A pa. 



pa 



-a(3(c 2 pp + C i 



(C-6) 
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A"^ = -a(5(c pp C pu + C%) , (C-7) 
A^ s = -a/3 (2C%) , (C-8) 



where 



[l-(3(S-q)Y 



Cpa — — — — — 2 j (C-9) 



c » = c " + i-KS- t ) - (C40) 

The elements of the submatrix B are expressed as 

B w = "« 2 /? 3 « (<t p k° p l) 2 - (a^) 4 » , (C-ll) 

= -a 2 (3 3 « (^a pL )(a pX ) 2 - (a pK f » , (C-12) 
#»vr* = -a 2 /? 3 « (a pX ) 4 - {<t p k) 4 »= , (C-13) 

where (• • •) denotes the average with respect to the weight appearing as the exponential function 
in (CT). This average should be evaluated at the replica-symmetric point to check the stability of 
the replica-symmetric solution against replica-symmetry breaking. The off-diagonal element of the 
block matrix (|C-2|) is given as 



S^ s = a/3 (SpySrt + SpsSay) . (C-14) 
We follow the usual procedure to find the eigenvalue of the replicon mode, 

Q P a = q + ri pa , (CT5) 
r pa = r + xrj p(7 , (C-16) 

where 

Vpv =V (P,<r + 1,2) , 
m = -(2-n)(3-n)r? , 

v PP = o • 

The eigenvalue equations are then given by 

J2 (A^ 8 + x5 pu ^ 5 ) i h/S = At? , (C-17) 
(xB pCT ~ iS + 5 pa ^ s ) r? 7 5 = Xxrj , (C-18) 
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with p, a 7^ 1, 2. Substitution of the expressions of A pa, ^ s , B pa, ^ s , d^' 1 ^ and into these equations 
leads to 

1 



x = A + 



[1-/9(5-9)1 



2 ' 



a/3 2 <C {ia pK <Tp L ) - (cj p ^) 2 ) 2 » +A 



1 , 



with A = A/a/?. This set of equations is solved for the eigenvalues as 



A± = ~^( u + v)± \J~^( U + v) 2 + l-uv 

where 

u = a(3 2 <C (((TpKO-pL) ~ {v p k) 2 ) > , 



" [l-P(S-q)} 2 q' 

Here we have used the equation of state (3.27) in the last equality. The AT line is determined by 
the vanishing point of the eigenvalue A+,& or equivalently, uv = 1, which has the following explicit 
form, 

a/3 2 « ((<J pK a pL ) - (a pK ) 2 ) 2 >= ^ . (C-19) 

It is necessary to evaluate the expectation values of spin variables appearing in ( CT9p . We define 
the quantity in the brackets <C • • • 3> in ( |C-19| ) as V: 

V = < ({a pK a pL ) - {<t p k) 2 ) 2 » 
= < (v p kVpl}(v\kv\l) > 

-2 < {<T p Kcr p L){cr\K){v*K) > 

+ < (<TpK) (c\k) (<7*k) {<?kk) > • 

According to the usual replica formalism,© this equation is equivalent to the following expression 
in the limit n — > 0: 



x cxp 



/3 x - a/3 2 \ - 

m-^}_^a pK + B}^a pK a PtK+1 + j^r ^ o pK a aL 



M 



Kp 



Kp 



KLpa 



» , 



(C-20) 



where different suffixes p, A, v and k correspond to different replicas. 

Let us consider the single-retrieval case m p = 5ni ■ m. The double summation over K,L,p,a in 
(C-20) can be decoupled using the Gaussian integral. In the limit n — > 0, we find 

, 2 



V 



D*(< 



<?kol) 



(C-21) 
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where (• • -) z stands for the average with respect to the weight 



exp 



1/ 



(C-22) 

K K K 

The averages appearing in the above equation can be calculated easily. The single-spin average is 
given by 

Tr a e Pha z +/3Aa x 



<j^. e /3hiT z +f3Aa x 
h 



t&nh(3\/h 2 + A 2 



V/i 2 + A 2 

with h = m + y/arz, and the first equality holds in the limit M 



(C-23) 



oo. The two-spin expectation 



value in (C-21) is the correlation function of the one-dimensional Ising model in a uniform magnetic 
field, and thus can be calculated by the transfer-matrix method. The result is, in the limit M — > oo, 



. h 2 A 2 cosh/3u(l - 2y) 

{vkvl)z = -o H 2 rTfl 

f z f z cosn pv 



(C-24) 



with w = (h 2 + A 2 ) 1 / 2 and y = (K - L)/M. 

The last expression ( C-24 ) of the correlation function has dependence on the Trotter numbers K 
and L through y, though we have calculated it under the static approximation which ignores such 



dependence. This inconsistency is remedied if we average ( C-24 ) over the interval < y < 1. For 
example, the ground-state expression of the order parameter S(= q) in the static approximation 
(3.33) can be obtained by the integral: 

,2 



lim 



Bz- 



(m + y/arzY 



(C-25) 



(m + \farz) + A 2 

This result seems quite natural in consideration of the definition (3.14) of the order parameter 
S(KL). Therefore, we replace the quantity in the brackets in ( |C-21| ) by its average, 

h 2 



J o dy {y K a L ) 



v 2 cosh 2 Pv ' Pv 3 
TA 2 

-•r (r-o) 



A 2 

+ —-it tanh Pv 



Equation ( |C • 21| ) now can be written as 



v = T 2 A 4 / 



and the AT line is finally obtained from this relation and ( |C-19| ) as 



q = arA 4 



Bz 



[m + yfarz) 2 + A 2 
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